Long-range correlation of neutrino in pion 

decay 

K. Ishikawa and Y. Tobita 
January 11, 2013 



Department of Physics, Faculty of Science, 
Hokkaido University Sapporo 060-0810, Japan 

Abstract 

Position-dependent property of the neutrino produced in high- 
energy pion decay is studied theoretically using a wave packet for- 
malism. It is found that the neutrino probability has a universal 
finite-distance correction if the pion has a long coherence length. This 
correction has origins in the interference of the neutrino that is caused 
by the small neutrino mass and a light-cone singularity of the pion and 
muon system. This term is positive definite and leads excess to the 
neutrino flux over the in-coherent value at macroscopic distance of 
near detector region. The flux decreases slowly with distance in a uni- 
versal manner that is determined by the absolute value of the neutrino 
mass and energy. Absolute value of the neutrino mass is expected to 
be measured indirectly from the neutrino detection probability at fi- 
nite distance. 
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1 Introduction 



Coherence phenomenon of the pion decay process becomes important when 
the pion is a wave of the long coherence length. The decay process then is 
studied using a time-dependent Schrodinger equation, and the neutrino in the 
final state becomes a wave that maintains the coherence for a long interval. 
Furthermore if this neutrino wave is a superposition of waves of particular 
complex phases, this shows an interference phenomenon. The total number 
of the neutrino events is large in the recent ground experiments, and the 
interference of the neutrino may become observable there. To study the 
interference phenomenon of the neutrino in the decay of high-energy pion 

ir ^ (x + Uf,, (1) 

especially the physics derived from the time-dependent wave function and 
sensitive to absolute value of the neutrino mass, is the subject of the present 
paper. 

Understanding the neutrino is in rapid progress. Neutrinos are very light 
and the mass squared differences were given from recent flavour oscillation 
experiments [H 121 [21 IH El E] using neutrinos from the sun, accelerators, re- 
actors, and atmosphere as, [7] 

Aml-^ ='m\-m\ = (7.59 ± 0.20) x 10"^ [eVVc% 
\A'm\^\ = \ml - mil = (2.43 ± 0.13) x 10"^ [eV^/c% (2) 

where rrii {i = 1 — 3) are mass values. The absolute values of masses are 
important but are not found from the oscillation experiments. Tritium beta 
decays [8] have been used for determining the absolute value but the exist- 
ing upper bound for the effective electron neutrino squared-mass is of order 
2 [eVVc^] and the mass is 0.2 — 2 [eV/c^] from cosmological observations [H]. 
The neutrino is also idealistic for the test of quantum mechanics due to its 
weak interaction with matters, although the validity of quantum mechanics 
is certain from many tests and applications made in the past. For these 
implications, much works should be done further. 

To study the coherence phenomenon of the neutrino, coherence properties 
of the pion should be known. One particle states which have finite coherence 
lengths are described by wave packets [101 [HI [12] • By using wave packets. 
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position-dependent amplitudes and probabilities are also computed, n Co- 
herence properties of the neutrino produced in pion decay are determined by 
the pion's coherence properties and decay dynamics. When the pion's co- 
herence length is short, the pion and the neutrino behave like particles. The 
total decay rate and average hfe time are computed in the usual way and the 
neutrino shows no physical phenomenon which has an origin in the spatial 
coherence. In other situation when the pion has a long coherence length, the 
pion and the neutrino behave like waves and the interference phenomenon 
of neutrino appears in the space-dependent wave function and probability, 
which will be presented here. 

We focus on quantitative analysis of the neutrino produced in pion decay. 
If these particles are described by plane waves, the transition amplitude has 
the delta function of the energy and momentum conservation and the total 
transition probability is proportional to time interval T. The decay rate and 
the mean life time are found from the proportional constant. Here we study 
finite-size corrections in the pion decay process. Namely the various prob- 
abilities of the pion decay process are computed in the situation where the 
distance between the initial and finial states is finite and the time interval be- 
tween them is also finite. From the integral representation of the probability, 
as is explicitly given in section 4, a large finite-size correction is lead if the 
correlation function of the pion and muon system has a light-cone singularity. 
The dependence of the probabilities on these values are also determined by 
the light-cone singularity. The time-dependent wave function of the whole 
process gives rise to this correction, which, furthermore, becomes observable 
quantity by the use of wave packet expression of the particle states, instead 
of the plane waves. 

A neutrino produced in pion decay propagates finite distance before it 
is detected. The distance is not fixed but varies within the size of the pion 
decay region. So the neutrino wave at the detector is a superposition of those 
that are produced at different positions. If these positions are in one pion, 
as in Fig. ([1]) the interference phenomenon of the neutrino wave is possible. 
Let the pion's coherence length, ^/a^, and velocity, v-^, and find the criteria 
when this happens. When the neutrino is produced either at time ti or t2 
from the pion that is prepared at Tt, and is detected at T^y, after traveling 

^The general arguments about the wave packet scattering are given in [T5 1 [T^ [T51 116) . 
In these works, however, situations where wave packet effects are neghgible were considered 
mainly. 
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with the hght velocity, they satisfy 



|(c(T, - ti) + v^{t, - T^)) - (c(T, - ts) + v^{t2 - T^))\ < V^, (3) 

if these particles travel in one- dimensional space, where the light velocity is 
used for the neutrino velocity Vi, = c. So this is one of the conditions for 
the neutrino interference in the one-dimensional space. If this criteria holds, 
the neutrino wave from ti can interfere with the neutrino wave from t2- At 
high energy, the pion velocity is close to the light velocity and the left hand 
side of Eq. becomes to c(m^/2i?^)(ti — ^2), hence this condition is sat- 
isfied at the time difference c{ti — ^2) ~ \f^{'^El / m^) . When this length 
c(ti — ^2) is macroscopic size, the coherence phenomenon of the neutrino at 
the macroscopic length is possible. We will estimate the coherence lengths 
of these particles in second section and will see that this condition in three- 
dimensional space is satisfied in the macroscopic distance. In this situation 
the amplitude for detecting neutrino at a finite distance is given by the inte- 
gral of the product of wave functions of the pion, muon, and neutrino over 
the coordinate. Now we recall that the velocity of any relativistic waves ap- 
proaches the light velocity, when the momentum becomes infinite. Hence the 
space-time positions of weak interaction i = 1,2 where the neutrino 

is produced in a system of the incoming pion wave and outgoing muon wave 
have the light velocity at the large momentum and the two point correla- 
tion function defined from the pion and muon states in the decay probability 
'-^■K,fi{ti — hyXi — X2) has a light-cone singularity. Since the light-cone sin- 
gularity is real and extended in wide area of space and time which almost 
overlaps with the neutrino's space-time path, the neutrino probability, which 
is the integral of the neutrino wave function multiplied by the light-cone sin- 
gularity, gets the direct effect from the phase of the neutrino wave function. 
The phase of the neutrino wave of mass mj, and energy is expressed by 
using the differences of two positions Ax and of two times At, as exp 
where is defined by EyAt — ■ Ax. This phase is expressed in the form 
(j) = {m1/2Ey) X cAt, if Ax = cUyAt is substituted. Here is the unit 
vector along the neutrino momentum. Since the angular velocity ml/2Ey 
is extremely small, this gives a large and long-distance correction. Thus the 
light-cone singularity of the pion and muon system and the light velocity and 
slow phase of the neutrino wave function are combined to produce the finite- 
distance correction to the neutrino probability. Because the neutrino's mass 
is extremely small, the correction is expected to be of the large distance. 
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Fig. 1: The geometry of the neutrino interference experiment. The neutrino 
is observed by the detector at and produced at t\ or ^2 in the inside of the 
travehng pion wave. 



We investigate the physical phenomena that have a deep connection with 
the neutrino's coherence and interferences in the high-energy pion decay. 
For simphcity we study the parameter region where the absolute value of the 
neutrino mass is much larger than the mass differences and the effect of the 
flavour oscillation is negligible. More general situations will be studied later. 
Other interference phenomena caused by solar neutrinos, reactor neutrinos, 
and others are studied in a next work. 

This paper is organized in the following manner. In section 2, the sizes 
and shapes of wave packets are studied. In section 3, we study neutrino 
production amplitude in the pion decay and in section 4, we study neutrinos 
detection probability. The probability is expressed as the sum of normal term 
and new term. The position-dependent probability is obtained from the light- 
cone singularity of the decay correlation function in section 5. Summary and 
prospects are given in section 6. 
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2 Wave packets 



In the decay of the pion Eq. ([T]), the pion is in the beam and the neutrino is 
detected by the detector. Pions are produced in colhsions of protons with nu- 
cleus and are superpositions of plane waves. Weak decay of a pion is studied 
normally using momentum eigenstates of the pion, neutrino, and muon. In 
the present paper, position dependent probability for the neutrino is studied 
and the neutrino is expressed by a wave packet, which is extended around its 
centers in position and momentum. In order to clarify this probability, initial 
pion is also expressed by a wave packet. The wave packet size is determined 
by the detection process for the neutrino but by the production process and 
medium effect for the pion. 

A particle in matter is described by a wave that maintains coherence 
in a finite region defined by the mean free path [11]. When the particle 
of the finite mean free path is emitted into vacuum, the wave representing 
this particle propagates. This wave is expressed by a wave packet. Hence a 
quantum mechanical transition of this particle state is computed from the 
overlap integral of the wave functions. The wave function of the finite spatial 
size has also a finite momentum width. 

We estimate the wave packet size of a proton first and those of a pion 
and neutrino next, following the method of our previous works (TUl ITT] . 

2.1 Size of pion wave packet 

Since pions are produced in collisions of protons with nucleus, the coherence 
lengths of pions are determined by those of a proton and nucleus. A proton 
in matter interacts with nucleus and has a finite mean free path or a finite 
coherence length. The target nucleus has a microscopic size of order 10~^^ [m] 
and its position is extended in a size of nucleus wave function in matter. Its 
magnitude is slightly larger than a nucleus intrinsic size. We use in the 
present paper the size of the order of 10~^^ [m] for the nucleus size. 

2.1.1 Proton mean free path 

The mean free path of a charged particle is determined by its scattering with 
atoms in matter by Coulomb interaction. An energy loss is also determined 
by the same cross section. Data on the energy loss are summarized well in 



6 



particle data summary [7] and are made use of for the evaluation of the mean 
free path of the proton. 

The proton's energy loss rate at the momentum, l[GeV/c], for several 
metals such as Pb, Fe, and others are 

dE 

- — = 1 - 2 [MeVg^^cm^l, (4) 
dx 

hence we have the mean free path of the 1 [GeV/c] proton in the material of 
the density p, 

E 1 [GeVl r n . ^ 

i^proton = ^ = J,, -1 , ^ = 50 - 100 cm . (5 

^ X p (1 — 2) X 10 [MeV g cm^g cm-'^J 

At a lower energy, 0.2 [GeV/c], the energy loss rate of the proton is about 
10 [MeVg^^cm^] and the mean free path is 

-^^proton = 10 [cm]. (6) 

The proton maintains coherence within the mean free path. Hence we 
use the mean free path for the coherence length of the proton feproton, 

(^•^proton -^proton • (7) 

When the proton of the finite coherence length is emitted into the vacuum 
or a system composed of dilute gas from the matter, it has the coherence 
length defined in matter. The proton has a constant coherence length in 
vacuum or in dilute gas when it is moving freely. The coherence varies when 
the proton is accelerated. If the potential energy V is added to the particle 
of momentum pbefore, then the momentum becomes Pafter and satisfies 



^Pbefore + + V = sj vl,,, + . (8) 

From Eq. (|H]) the variants of the momentum satisfy 

■^efore ^ "^Pbefore = "^after ^ alter i 

Pbefore Pafter /q\ 

■^^before = , „ „ i "^aftcr — ; r, ^ • v"/ 

V Pbefore + ^ V Pafter + ^ 

Hence the coherence length of a particle, (Jxbefore, which is proportional to 
the inverse of ^Pbefore^ becomes to feafter after the acceleration from a velocity 
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Vheforc to cL velocity f after- The Coherence length is determined by the velocity 
ratio, 



(^a^after = ^a^before X (IQ) 
^before 

The velocity is bounded by the light velocity c, and the velocity ratio from 
1 [GeV/c] to 10 [GeV/c] is about 1.2 and that from 0.2 [GeV/c] to 10 [GeV/c] 
is about 5. Hence the proton of 10 [GeV/c] regardless of the energy in matter 
has the mean free path 

^Xproton ^ 40 - 100 [cm], (11) 
in vacuum or a dilute gas. 



2.1.2 Pion mean free path 

Coherence length of pions which are produced by a proton collision with 
target nucleus is determined by the proton's initial coherence length Eq. dlip 
and the target size 10~^^ [m], which is negligibly small. Since the time interval 
of the proton wave is the same as that of the pion wave, the coherence length 
of the pion, 5xpion, is given from that of the proton, feproton, in the form 

, fu^pion "■^proton ~ "■^proton- \ J 

^proton ^pion ^proton 

In relativistic energy region, particles have light velocity. Consequently from 
Eq. (fTTI) . we have the pion's coherence of 1 [GeV/c] or larger momentum 

Sxpion - 40 - 100 [cm]. (13) 

We use this value of Eq. ( [T3l) as the size of the wave packet 

y/On = SXpion, (14) 

in latter sections. 

In vacuum and dilute gas, pions propagate freely with the same coherence 
lengths. From Eqs. ([5]), ( ITTl) . (|T3l) . the proton and pion have the coherence 
lengths of the order 50 — 100 [cm]. 
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2.2 Size of neutrino wave packet 

A size of wave packet for observed neutrino is determined from its detection 
process. Neutrinos interact with nucleus or electrons in atoms incoherently 
and the neutrino event is identified by the reaction process in detector. Hence 
the size of the observed neutrino is not determined by the mean free path 
but by a size of the minimum physical system that the neutrino interacts and 
gives the signal. They are either nucleus or electrons in atom. The nucleus 
have sizes of order 10~^^ [m] and the electron's wave functions have sizes of 
order 10~^° [m]. So neutrino wave packet is either 10~^° [m] or 10^^^ [m]. 
The interactions of the muon neutrino in detectors are 

z^/. + + i/^, (15) 

'^/i + /i" + z/e, (16) 

+ + {A + l) + X, (17) 

+ u^, + A + X, (18) 

hence the size of the neutrino wave packet in processes ( ITSi) and ( fT6|) is 
of order 10~^°,[m] 

-10 



^=10-^"[m], (19) 



and the neutrino wave packet y/d\, in processes ( lT7j) and (jT8|) is of order 

10-15 [m] 

^=10-iMm]. (20) 

The interactions of the electron neutrinos in detectors are 

i/g + — )■ + z/e, (21) 
Ve + A^e- + {A + l)+X, (22) 
Ve + A^ e + A + X. (23) 

The neutrino wave packet ^ov in processes (|2T]) is of order 10"^° [m], Eq. ( IT9l) . 
and the neutrino wave packet -y/ov in processes f l22|) and fl23|) is of order 
10~i5 [m], Eq. (120 p . They are treated in the same way as the neutrino from 
the pion decay. 

From Eqs. (I19p and fl20|) . the neutrino have the wave packet sizes of the 
order 10"^° [m] or 10"^^ [m]. 
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We study neutrinos described by the wave packets of these sizes in many 
particle processes. In this respect, the neutrino wave packet of the present 
work is different from some previous works of wave packets that are connected 
with flavour neutrino oscillations [171 UHl |T9l [201 |2T1 [221 123] , where one particle 
properties of neutrino at production are studied. It is important to study 
the neutrino wave packet for our purpose of studying the interference. 

2.3 Wave packet shape 

A particle of the finite coherence length is described by a wave packet, which 
is extended in the momentum and position around the centers {p,X). Al- 
though its precise shape is unknown generally in real experiments, the physi- 
cal quantity should have universal value that is independent from the details 
of wave packet. We shall obtain the physical quantity that has universal 
properties. The quantity that depends on the wave packet shape is neither 
genuine nor universal and is not important. 

We give a summary of the wave packets in the following, since the wave 
packet may be unfamiliar to some readers. First, the wave packets are lo- 
calized in the momentum and position around their centers p!3| [TH [151 [16] . 
The whole set of the wave packets becomes complete set [TO] . 



independent from the shape parameter /3 which is discussed later. 

Second, the wave packet preserves the discrete symmetries such as invari- 
ances under space and time inversions, which have not been studied before 
but are quite natural, since the origin of the wave packet is the interactions 
of the particle with matters in detectors and the wave packet sizes are esti- 
mated in the next section under this consideration. This interaction has an 
origin in quantum electrodynamics that preserves parity and time reversal 
symmetries. So the wave packets should preserve parity and time reversal 
invariances. We study, hence, the wave packets which are superpositions of 
the plane waves around the central momentum with a weight function that 
has the same property under these transformations. 




(24) 




(25) 
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where the momentum p is the central value of the momentum. In the present 
work, k is used for the integration variable and p is used for the central value 
of momentum. 

Since the time reversal invariance is satisfied in the quantum electrody- 
namics and weak interactions in the lowest order, the wave packet which 
is invariant under the time inversion is most important. Under the time 
inversion, the coordinate and momentum variables are transformed into 

X — >■ X, 

— * — * 

k — >■ —k, p — >■ —p 

t -t, (26) 
and the plane wave is transformed to its complex conjugate, 

^i{Et—k-x) ^ ^—i{Et—k-x) ^^i{Et—k-x)^* (27) 

So when the weight satisfies 

w{k;p) ^ {w{-k; -p)) = {w{k-p)y, (28) 

the state described by the wave packet is the superposition of the waves 
which are transformed equivalently under the time inversion, 

j dkw{k]p)e'^^^-^-^^ ^ J dkw{-k;-p){e'^^'-''-'^^)* 

dkw{k;p)e'^^*-^-^A . (29) 



Next we study the space inversion. Although neutrino violates this sym- 
metry, the wave packet may preserve. Under the space inversion, the coor- 
dinate and momentum variables are changed into 

— * — * 

X — >■ —X, k — >■ —k, p — >■ — p, (30) 
and the plane wave is changed into 

^i{Et-k-x) _^ ^i{Et-k-x)_ ^gj^^ 

So when the weight satisfies 

w{k;p) ^ w{—k; —p) = w{k;p), (32) 
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the state described by the wave packet is transformed in the following way, 
j dkw{k-p)e'^^'-^ '^^ j dkw{-k--p)e'^'^'-^-'^^ 

= j dkw{k]p)e'^^'-^-^\ (33) 

in the same way as the plane wave. 

The invariance under the time reversal, Eq. (!28|) . is a strong condition 
and leads the important result. We study wave packets that satisfy Eq. ( !28|) . 
The simplest form of satisfying this property is the Gaussian wave packet 

|p-X,/3o) = ^ / rf^e-f(^^-^^)'e<^(^^)(*-^)-^^-(^^-^)), (34) 
(27r)2 J 

where the parameter a shows the size of the wave packet in the coordinate 
space and is the normalization factor. We write a wave packet of a pa- 
rameter [3 as |p, The Gaussian wave packet is the state at /3 = /3o 
and non-Gaussian wave packets are the states at /3 7^ /3o- The calculations 
using the Gaussian wave packet are presented in most places hereafter for 
the sake of simplicity, except the derivation of the large distance behavior of 
the neutrino wave function. It has the universal form regardless of the details 
of the wave packets, as far as the invariance under the time inversion holds. 
The theorem for the general wave packets that are invariant under the time 
inversion Eq. fl28|) holds and is given in Section 4. It is worthwhile to verify 
the results using explicit calculations of the Gaussian wave packet or various 
non-Gaussian wave packets, which are made partly in the text and in the 
appendix. 

The normal physical quantity of microscopic physics is obtained from 
ordinary scattering which has no dependence upon the distance or time in- 
terval between the initial and final states. This is because the length of the 
microscopic system is so small that the size of experimental apparatus is 
regarded infinite and the boundary conditions of ordinary scatterings which 
are defined at t = ±00 of plane waves are suitable. Hence the amplitude 
and probability are defined using plane waves and the boundary conditions 
at t = ±cx3 ensure the independence of the probability from the distance and 
particle's coherence length. 

For the neutrino the situation is different because the neutrino mass is 
so small that a new energy scale defined by irt^/Ey becomes extremely small 
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and a spatial length which is inversely proportional to this energy becomes 
macroscopic length. It is important to find a physical quantity that depends 
on this length, since it may give a new important information. To find a 
scattering amplitude and probability that have a dependence on the length 
or time interval, we use the amplitude defined from wave packets of finite 
spatial sizes, a. Those wave packets that have central positions and are 
localized well around the central positions and that have the same proper- 
ties in the momentum variable are suitable for this purpose. The simplest 
wave packets of having these properties are Gaussian wave packets Eq. (IMI) 
which satisfy the minimum uncertainty relation between the variances of 
the position and momentum. Non-Gaussian wave packet parameterized by 
a parameter /3 and has larger uncertainties is easily defined by multiplying 
Hermitian polynomials to the Gaussian function. 

From Eq. fl24|) . the total transition probability is independent from the 
parameter (3 

I {neutrino; p^, X^, (3; muon\T\pion; pV, ^tt) |^ 
= I (neutrmo; Xy, /3o; w^^^o^^|T|p^or^,; Ptt, -^715 ^7r)p5 (35) 

which also agree to that of the plane waves. The probability of the finite 
distance is defined by restricting the center positions X in the inside of a 
finite spatial region V . The completeness condition for the function in this 
region V is then 

|p,X,/3o)(p,X,/5o|= \P.X,P){p,X,Pl (36) 

p,x<y p,x<v 

and the probability satisfies 

I {neutrino; py, Xy, (3; muon\T\pion; pV, ^n, ^7r)P 

Pv,Xv<V 

= Yj \{neutrino;pu, X^,, (3Q;muon\T\pion;p-„, Xt,,T.„)\'^ , (37) 

Pu,X^<V 

and is also independent from the (3. We study the probability that is the 
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average over a finite neutrino energy region Vp 

\{neutrino;p,,,X,y, (3;muon\T\pion;pT,,XT,,TT,)\'^, (38) 

p„<Vp,X„<V 

and extract the physical quantity. For the genuine physical quantity, the 
value should be determined uniquely. The energy uncertainty of the neu- 
trino experiment is of the order 10 percent of the total neutrino energy. This 
magnitude is the same order as that of the minimum uncertainty if the neu- 
trino energy is 1 [GeV]. So the minimum wave packet is applied directly 
in this case. The probability for a larger energy uncertainty is computed 
using the probability of the small energy uncertainty. The dependence of 
the neutrino probability on the distance will be shown to have the universal 
behavior. 



3 Position-dependent amplitude of neutrino 

Using the wave packet representation obtained in the previous section, we 
find the position-dependent neutrino amplitude. 



3.1 Semileptonic decay of the pion 
3.1.1 Weak interaction 

Semileptonic decay of a pion is described by the weak Hamiltonian 

Hw = g J dxd^Lp{x)Jy_^{x) = -igrrif^ J dx(f{x)J5{x), 
J^-Ai^) = /i(x)7^(l - l^Hx), Mx) = /i(x)(l - ^^Hx), (39) 

where (p{x), n{x), and z/(a;) are the pion field, muon field, and neutrino 
field 121]. In the above equations, g is the coupling strength, Jy_y^{x), and 
J5{x) are the leptonic charged vector current, and leptonic pseudoscalar. The 
coupling is expressed with Fermi coupling Gp and pion decay constant f^, 

9-%U. (40) 
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3.1.2 Neutrino production amplitude in the pion decay 



When the number of the neutrino events is large enough, the distributions 
of the neutrino events are computed with the quantum mechanical wave 
functions. For the electron bi-prism experiments of Tonomura et al |25] . 
the electrons' interferences of the sharp energy become visible when the to- 
tal number becomes significant. Even though the initial electrons have no 
correlation and each event occurs randomly, when number of the electrons 
becomes large, the signal exceeds the statistical fluctuation and the distri- 
bution shows the interference behavior. A large number of events is also 
necessary for the quantum mechanical interferences become visible for the 
neutrinos from the pion decays. In this situation each pion decay occurs 
randomly and its distribution is computed using quantum mechanical wave 
functions. We study this situation hereafter and compute the amplitude and 
probability quantum mechanically. 

A wave function which describes a pion and its decay products satisfies a 
Schrodinger equation 



where Hq stands for the free Hamiltonian and H^, stands for the interaction 
Hamiltonian Eq. (!39|) . The solution is 



in the first order of H^. A muon and a neutrino are described by one wave 
function of Eq. (^2\i . |muon, neutrino (t)), and all the informations of this sys- 
tem are computed from this wave function. Physical quantities were studied 
with this wave function [23] and the finite time interval effect was found. 

Here the matrix elements between |muon, neutrino (if:)) and the state u) 
defined next is studied instead of the wave function itself. Since all the infor- 
mations of the wave function |muon, neutrino (t)) are included in this matrix 
element, as far as a complete set of \fi, v) is used, this matrix element is 
equivalent to the wave function. Furthermore the position-dependent infor- 
mations on the decay process are found by the use of the wave packets. The 
coherent phenomenon of the neutrino which depends on the position and the 
neutrino detection probability are studied here with this matrix element. 




(41) 




(42) 



15 



The transition amplitude from a pion located around the position X^^ to 
a neutrino at the position and the muon is studied hereafter. Because the 
muon is not observed and all the states are summed over, the most convenient 
wave is used for the muon. This amplitude is 

T = j d'^x{ii,u\H^{x)\7r), (43) 

where the states are described by the wave packets of central values of mo- 
menta and coordinates and the widths in the form 

\tt) = \p^,X^,T^), = \fi,p^;u,pi,,Xi,,T^). (44) 

The states are defined by the matrix elements, 

3 

z/,p^,X^,T^|/i(x)(l - 75)z/(a;)|0) 



1 



Ny f .r '^i^ it ^ \2 I m,, \ ( m 



(2vr)3 7 \Eip,)J \Eih 



where 



3 , X 3 / I 



In this paper the spinor's normalization is 

E 'j ■ p + m 
u{p, s)u{p, s) = — — — . (48) 

s 

In the above equation the pion's life time is ignored. The sizes, 0",^ and a^, 
in Eqs. (1451) and (H6|l are sizes of the pion wave packet and of the neutrino 
wave packet. Minimum wave packets are used in most of the present paper 
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but non-minimum wave packets are studied and it is shown that our results 
are the same. 

From the results of the previous section, the size of pion wave packet 
is of the order 0.5 — 1.0 [m] and the momentum has a small width. So the 
pion momentum is integrated easily, and is replaced with its central value 
and the final expression of Eq. fHSl) is obtained. For the neutrino, the size of 
wave packet should be the size of minimum physical system that a neutrino 
interacts, i.e., the nucleus. Hence to study neutrino interferences, we use the 
nuclear size for cxy. 

The amplitude T for one pion to decay into a neutrino and a muon is 
written in the form 

T = J dtdxT{t,x),T{t,x) = j dKT{t,x,K) (49) 

where 



a J \E,EJ (27r)3 

This amplitude depends upon the coordinates (t, x) explicitly and is not 
invariant under the translation. 



3.2 Neutrino wave function 

We integrate the neutrino momentum of Eq. by applying Gaussian in- 
tegral and find the coordinate representation of the neutrino wave function. 
It is found that a phase of neutrino wave function has a particular form that 
is proportional to the square of the mass and inversely proportional to the 
neutrino energy. 

^For the non-minimal wave packets which have larger uncertainties Hermitian polyno- 
mials of ki, — -pv are multiplied to the right-hand side of Eq. (|46|) . The completeness of 
the wave packet states is also satisfied for the non-minimum case and the total probability 
and the probability of the finite distance and time are the same. We will confirm in the 
text and appendix that the universal long-range correlation of the present work is inde- 
pendent from the wave packet shape as far as the wave packet is invariant under the time 
inversions. 
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3.2.1 Phase of neutrino wave function 



For not so large t — T^ region, ky is integrated around tlie central momentum 
Pu in Eq. poj) . The amplitude becomes, 



where N is the normalization factor, Vy is the neutrino velocity, and is the 
phase of neutrino wave function. They are given by 

3 3 

27r\ 2 /27^^ ^ 



'^^'^-'^"Uj '''' 

= E{p,){t - n) - p„.(f - .?„). (63) 

The phase 4>{x) is rewritten for a small wave packet by substituting the 
central value x of neutrino's Gaussian function 



x = X, + Vy{t-n), (54) 



in the form 



which has a typical form of the relativistic particle. The phase becomes 
proportional to the neutrino mass squared and inversely proportional to the 
neutrino energy. The magnitude of the phase at the position Eq. ( 15^ is 
small due to the cancellation between the oscillation in time and space. This 
cancellation does not occur in the derivatives of the phase with respect to 
the coordinate, which are given in the form 

^ (jy = p^y. (56) 



dXf, 



This is not proportional to the square of neutrino mass but is determined by 
the energy and momentum. 

When the position is moving with the light velocity 

x = X, + c{t-T,), \c\ = 1, (57) 
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then the phase is given by 



2 

mf. 



and becomes a half of (pg. 

Interference of the neutrino due to this slow phase of the energy E^, and 
El, + AEi, is identical if the neutrino phase 0c satisfies 

d - 

A(f), = 5E,-—(j), < TT. (59) 

Hence the neutrinos of the energy width AE^ show the same interference. 
We will see later that this is actually satisfied for the light neutrino of the 
energy width around 100 [MeV]. 

The slow phase 0c of Eq- dSS]) is a characteristic feature of the neutrino, 
and is the intrinsic property of the neutrino wave function along the light 
cone. Using the Gaussian wave packet, this phase has been obtained easily. 
Since this phase is genuine of the neutrino wave function, the same behavior 
is obtained in more general cases including the spreading wave packet or 
non-Gaussian wave packets. The spreading of the wave packet is negligible 
in the longitudinal direction and is not so small in the transverse direction. 
Its effect has been ignored for simplicity in this section and will be studied 
in the latter section and appendix. It will be shown there that the spreading 
of the wave packet in the transverse direction modifies the k^, integration but 
the final result turns actually into the same. For general wave packets, if 
the completeness, Eq. (l36l) is applied, the total probabihty satisfies Eq. ( 137|1 . 
The transition amplitude in this case has the same phase 0(x), Eq. (15 3 p and 
is written in the form, 

X n(p^)(l-75)zy(p;)e^<^(--)ti;(x-X,-i;,(t-T,)), (60) 

where w{x) is the amplitude of neutrino wave function obtained by the 
Fourier transformation from the w{p,k). The spreading effect is negligible 
for the small time interval and w{x) is given in the form 

w{x)= [ dke^^-^w{k), k = k-p, (61) 
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and decreases fast with Moreover from Eq. fl28p . this function is real 
function 

w{x) = w{xy, (62) 

which leads the universal long-distance behavior to the neutrino probabil- 
ity. The spreading effect is included when the long-range component of the 
correlation is studied in Section 4. 



3.2.2 Position-dependence, energy momentum non-conservation, 
and interferences 

When the space-time coordinates {t, x) are integrated in the amplitude of the 
plane waves, the delta function of the energy and momentum conservation 
emerges. The scattering amplitude with this delta function shows that the 
final states have the same energy and momentum with the initial state. On 
the other hand, the posit ion- dependent amplitude T for the wave packets is 
not invariant under the translation and has no delta function. So the energy 
and momentum of the final state are not necessary the same as the initial 
state. The states which do not satisfy the energy and momentum conser- 
vation should be included to get consistent results from the completeness. 
This amplitude shows the position-dependent behavior, from which a new 
information is found. 

Let us compare the neutrino velocity with the light velocity. The neutrino 
of energy 1 [GeV] and the mass 1 [eV/c^] has a velocity 

vlc=l- 2e, e = ' = 5 x 10-^^ (63) 

hence the neutrino propagates the distance Z, where 

/ = /o(l _ e) = /g - 5/, 6l = loX e, (64) 

while the light propagates the distance Iq- This difference of distance, 61, 
becomes 

5l = 5x 10-^^ [m]; Iq = 100 [m], (65) 
5l = 5x IQ-^^ [m]; Iq = 1000 [m], (66) 



which are much smaller than the sizes of the above wave packets Eq. f[T3|) . 
This gives the important effect for the neutrino amplitude at the nuclear 
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target or the atom target to show interference. The geometry of the neu- 
trino interference is shown in Fig. [T] The neutrino wave produced at a time 
ti arrives to one nucleus or atom in the detector and is added to the wave 
produced at t2 and arrives to the same nucleus or atom same time. A con- 
structive interference of waves is shown in the text. 

4 Position-dependent probability and inter- 
ference 

The probability of detecting neutrino at a finite distance is studied in this 
section. Particularly the finite-time interval T correction, i.e., the deviation 
of the transition probability from the T-linear form, is obtained. For this 
purpose, the total probability is expressed with the two point correlation 
function of decay process, which has the light-cone singularity. The light- 
cone singularity leads the finite-distance correction of universal form. 

4.1 Transition probability 

The transition probability of a pion of momentum located at space- 
time position {Tt^,Xt^) to a neutrino of momentum p,^ at space-time position 
(Tj,, Xj,) and a muon of momentum is written in the form 



where 5'5(si,S2) stands for the products of Dirac spinors and their complex 
conjugates. 




(67) 



S5{Sl, S2) = {u{Pf,){l - l5)v{Pv)) {u{p^){l - l^)v{pu)) 

and its spin summation is given by 



(68) 




(69) 
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For the plane waves, the coordinates Xi and X2 are integrated easily and 
{5^^\pf — Pi)y = S^^\pf — Pi)S'''^KO) is obtained. Thus the four dimensional 
momentum is conserved and the probability is proportional to the volume V, 
which is cancelled with the normalization factor and T from 5^^\0) = VT. 
The probability is proportional to T and the proportional constant gives the 
decay rate per time. Here we calculate the finite-time interval correction 
for the neutrino detection process. The neutrino is observed through the 
interaction with the nucleus and is expressed with the wave packet. In order 
to obtain the finite-time interval correction, we write the probability with the 
correlation function of the pion and muon system. The transition probability 
which is integrated in momenta of the final state and average over the initial 
momentum is written, then, in the form 



dpTrPexpiPn)dp^dpu ^ |T| 



2 



^ {2Tcy J 

/47r\ ^ /47r\ ^ / 1 V 

using the correlation function ATr^^{6t,6x). The correlation function is de- 
fined with a pion's momentum distribution PexpiPir), by 

{27vy J E{p^) E{p^) 

(71) 

In the above equation, Eqs. fl55l) and fl56l) were substituted. The muon 
momentum is integrated in whole positive energy region, because the muon 
is not observed. If the muon is observed and its momentum is measured 
together with the neutrino, then the muon momentum is integrated in the 
finite-energy region as in the inclusive hadron reactions [26] . Pion momentum 
is integrated in the region specified by the initial pion beam. The velocity {7^ 
in the pion Gaussian factor was replaced with its average v-j^. This is verified 
from the large spatial size of the pion wave packet discussed in the previous 
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section. If the pion's momentum is fixed to one value, the correlation function 

(72) 

is used instead of Eq. flTT]) . 



4.2 Light-cone singularity 

The correlation function AT^^^{6t,6x) has a singularity near the light-cone 
region 

X = 5t^- \Sxf = 0, (73) 

which is extended into a large \6x\ and is independent fromp^. So A^_^(5)f:, (5x) 
has the same light-cone singularity. So the probability Eq. ( |70l) has a large 
finite-distance correction due to the singularity of AT,^^{6t,6x). We find the 
light-cone singularity of A^.^(5t, 6x) [27] in the following and apply the result 
also to A^^^(5t, 5i;). 

4.2.1 Separation of singularity 

If the particles are plane waves, the energy and momentum are strictly con- 
served and the momenta satisfy 

Ptt = + P'^' (Pvr - P^lf = ml^O. (74) 

Hence the momentum difference p,r — P/i is almost on the light cone and the 
A.„^^{6t,6x) around the light cone, A = 0, is important in Eq.(170l). In order 
to extract the singular term from ATr^n{6t,6x), we write the integral in the 
form 

A.„(5Mx) = ^-^/(p.,5x), 



Pti-^t~^^ 



^-ii{E{p^)-E(pf,)}5t-{p^-p^)-Sx) 



(75) 
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Next the integration variable is changed from to q — — p^r that is 
conjugate to Sx. We have then 



I{p^,dx) = ^ j d'^(ld{q^ + pl){{p^ + q)-py)lm 



TT 



.(? + P7r)^ -ml-ie 
1 



r,iq-Sx 



+ d^qd{q''+pl){q-p,)\m 



{q + Pn)^ - ml - ie 
1 



iq-Sx 



.{q + PnY - ml - ie 
and we separate the integration region into two parts: 



JqSx 



(76) 



/(p^, 5x) = /i(p^, (5a;) + hiPn, Sx), 

d 
d6x 
1 



h{pTV: Sx) -- 



P-K-Pv+Pv- -I 



(q + Pn)'^ -ml-i 



le 



Jq-Sx 



hiPn^Sx) = - 

TT 



dq^(rqp,,-{pT, + q)lm. 



{q +P7r)^ - ml - ie 



^iq-Sx 



(77) 



Ii{Ptt^Sx) is the integral of the region Q < q^ and has a singularity and 
h{PTT,Sx) is the integral of the region — p° < Q'" < and is regular. Al- 
though the large finite-distance correction is derived from the singular term 
h{Pn, Sx), Ii does not contribute to the total probability at an infinite time 
for the plane waves. Conversely hipTv, Sx) contributes to the total proba- 
bility without finite-distance correction for the plane waves. Consequently 
by expressing / into the sum of Ii and I2, the correction at finite distance 
is found easily. So the physical quantity at the finite distance that refiects 
interference is computed using the most singular term of Ii. 
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4.2.2 Correlation function 



The denominator of the integrand of Ji is expanded in p^-q and the Ji is 
written in the form 



IT 



d^q9{q'^) Im 



1 



g2 _|_ ^2 _ ^2 _|_ 2q-p^ — ie 



d^qe{q^) Im 



q2 + m-^ - It 



2p^-q 



- I d\9iq')\l-2p^- 



d \ d 



d 



d 



d6x J dm? 



+ 



q2 _|_ ^2 _ 
Im 



^iq-Sx 

+ ■ ■ 

1 



iq-5x 



^ig-Sx 



= 2 <j 1 - 
where 



d6x J dm? 



+ 



d\e{q'')5{q' + m')e 



gZ _|_ _ 

~ 2x jg-fe 



~ 2 2 2 

m = m^ — m^. 



The expansion in 2q-p of Eq. ( 178|) converges in the region 

2pn-q 



g2 _j_ ^2 



< 1. 



(78) 



(79) 



(80) 



Here q is the integration variable and varies. So we evaluate the series after 
the integral and find the condition for the convergence. We find later that 
the series after the momentum integration converges in the region 



2pn-Pu 

m? 



< 1. 



(81) 



So the following result that is obtained using this expansion is applied in 
this kinematical region. In the outside of this region, the evaluation of the 
integrals Ji and I2 separately is not useful and I is integrated directly. 
The formula for a relativistic field of the imaginary mass 



d^qe{q^)5{q' +m')e 
im 



{2-Kfi 



An 



6{X)e{6t) + f short 



short 



im 



e{-X) |a^i (mV^^ - te{6t)Ji (mV^^ | 



(82) 
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Fig. 2: The region I corresponds to short-range correlation where t ~ x ~ 0. 
On the other hand the region II corresponds to long-range correlation where 
— ~ and both t and x can be macroscopic. 



where Ni, Ji, and Ki are Bessel functions, is substituted to Eq. (!78|) . The 
first term in the right hand side of Eq. (|82i) is the most singular term and 
the second and third terms have singularity of the form 1/A around A = 

and decrease as e~™"V^' or oscillates as e*™V^'. The singular functions and 
regular functions behave differently and are expressed in Fig. [2] for one space 
dimension. The singular functions have the value around the light cone and 
the regular functions have finite value in small area around the origin. Since 
the light cone is extended in macroscopic area, the light-cone singularity 
leads the correlation function to become long-range. The long-range correla- 
tion function from the light-cone singularity and the short-range correlation 
function from the regular function are computed next. Thus the correlation 
function Ji becomes long-range only along the light-cone region and decreases 
exponentially or oscillates rapidly in other directions. So then is 

written in the form 



dSx J div? I \ 47r 



short 



(83) 
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Next I2 is evaluated. For I2 we use the momentum q = q + p-^ and to 
write in the form 



d q (pj,-g)Im 



0<g'J<pO 



q2 _ jyfl _ 



d 



dSx 



^ ■/0<(jO<pO 

d?q 



m. 



d[pl-Jq^ + ml)e 



2 \ Ji-^x 



(84) 



The regular part I2 has no singularity because the integration domain is finite 
and becomes short-range. 

Thus the first term in Ji gives a finite long-distance correlation and the 
rests, the second term in Ji and J2, give short distance correlations. The 
correlation function, A^^^(5t, 5x) has a singular term and a regular term and 
is written in the form 



1 



1 



. d \ d 



Ptt-Pu+Pp 



. d 

' dSx 



2(27r)^i 



—6{X)e{6t) + f short ) + I2 



d6x J dm? 

where the dots stand for the higher order terms. 

4.3 Integration of spatial coordinates 

Next, the coordinates Xi and X2 are integrated in 

X A^,^(5t,(5x). 



^5) 



(86) 



The derivative in the above integral is computed using the integration 
by part as 

dx{p^ - A;,)e'(^(^")-P--^"V(a;), (87) 
where a function f{x) is an arbitrary function and Eq. (1561) was used. 
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4.3.1 Singular terms: long-range correlation 

The most singular term in Eq. (|86|1 is 



X ^5(A)e(5t) (88) 



and is rewritten using the center coordinate = ^ ^ ^ and the relative 
coordinate f = xi — X2 in the form, 



X ^6{\)e{6t). (89) 

The center coordinate X is integrated easily and Js^x) becomes the integral of 
the transverse and longitudinal component (ty, r;) of the relative coordinates, 

Finally this is computed in the form 



/ ^s.a^ 1 isjst) — 

- K^)^y^<5t)e^^^(^*). (91) 
The next term in Eq. (1861) is from 1/A. We have 
Ji/A = f di,df2e'*<'"e-5t(*'-*'-"-"'-^">)'e-jfc("-*'-'"<"-^'-')' 



which becomes 



1 



^"^ ''"^^^2{^J e--^e*-. (93) 

28 



This term has the universal \St\ dependence but its magnitude is much smaller 
than that of Js(x) and is negligible in the present decay mode. 

From Eqs. fl9T|) and f l93l) . the singular terms Js[\) and Ji/\ have the slow 
phase 4>c{St) and the magnitudes that are inversely proportional to the time 
difference. Thus these terms are long-range with the small angular velocity 
and are insensitive to the m^. These properties of the time dependent cor- 
relation functions J5(a) are most important for the neutrino probability to 
have the universal long-distance behavior and also hold for the general wave 
packets of Eq. (!60!) which satisfy Eq. (162|) . 

Theorem 

The singular part Js(_x) of the correlation function has the slow phase that 
is determined by the absolute value of the neutrino mass and the magnitude 
that is inversely proportional to the time difference at the large distance, of 
the form Eq. flOTl) . The phase is given by the sum of 4>c{St) and the small 
correction, which is inversely proportional to the neutrino energy in general 
systems and becomes if the neutrino wave packet is invariant under the 
time inversion and is real. 

(Proof: General cases including spreading of wave packet ) 
We prove the theorem for Js(x) which is written in the form, 

^^(A) = j dfe'^^^'^^wiir-vM)) X ^6{X)e{6t), (94) 

where w{x—vt) is the wave packet in the coordinate representation. Since the 
light-cone singularity is extended in large time, the spreading effect becomes 
important and is included in the time-dependent correlation function Js{^)- 
The wave function w{x — vt) with the spreading effect is expressed in the 
following form 



Cij — Co6ij, Co — (96) 



instead of Eq.f l6T|) . The quadratic term of k in the expansion of E yp+ kj is 

included and this makes the wave packet spread with time. The coefficient 
Cij in the longitudinal direction is negligible for the neutrino and is neglected. 
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w{ri — Vi,5t, tt) is written with this w{x — vt) as 
w{ri - vjt, vt) = J dXw {x +''-^w*[x-'- 

^ J dkidkTe''''^''-''"^^^+'''^-^^+'''°^^^^^ (97) 
Then the correlation function becomes 

J 47r( ^n! V< 

X 5{5t'' - rf)ei6t) 

dridrTdkidkTe''^^^''''^^''''^''-''-'^'^^^^^ 



J5{\) 

+2 2 



a(sty 



47r I ^ 

\ n=l 



X 7^ <! 1 + > : iii-r^r (^) } - rf)em 

dridkie'^^^'''''^+'^' -^•'^'^ drTdkre+'^'^'T^' | w ( A;; , A^r) T 



X - ^ 1 + E h^^ririT (^(Ij,)" e--'>^(«= - ^ym. m 



The variables tt are integrated first and kr are integrated next. Then we 
have the expression 



5(A) 



>^ s{^ + E;^(-«ft«r(^)"}(2W--aeW. (99) 



Using the following identity 



and taking the leading term in 1/St, we have the final expression of the 
correlation function at the long-distance region 



cm 
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Hence Js{x) in Eq. fllUip becomes almost the same form as Eq. ( IHTI) and 
the slow phase (pd^t) is modified slightly and the magnitude that is inversely 
proportional to the time difference. Js(^x) has the universal form for the 
general wave packets. By expanding the exponential factor and taking the 
quadratic term of the exponent, the above integral is written in the form 



/ 



dki{l + ki{i{l - v,)5t + Co) + ^{ki{t{l - v,)5t + Co)f)\w{ki, 0)|2 
wo{l + - v,)6t + Co) + - v,)6t + Co)') 



= wo(^ + Cod, + - (1 - v,f6t'^ + i{di{l - v,)5t + ^2^^0(1 - v,)5t). 

(102) 

We substitute this expression into the correlation function, and we have 

i^,(<5t)(l+<5) r 

J5(A) =7re-^''^u;o(l + 7)e(5t) , = / dki\w{kiM\ (103) 

where the correction terms are given by 

di d2 I di d2 [ I \^ , -2c.2 

d, = — [ dkiki\w{k,0)\^, d2 = — [ dkikf\w{k,0)\^. (104) 

Wo J Wo J 

In the wave packets of time reversal invariance, \w{ki, 0)^ is the even function 
of ki. Hence di vanishes 

di = 0, (105) 

and the correction are 
Q.E.D. 

The light-cone singularity is along — 5x'^ = 0, which is so close to the 



neutrino path that it gives a finite contribution to the integral Eq. fl86|) . Since 
the light-cone singularity is real, the integral is sensitive only to the neutrino 
phase and shows interference of the neutrino. 
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4.3.2 Regulcir terms: short-range correlation 

Next we study the regular terms. The regular term is short-range and the 
spreading effect is ignored and the Gaussian wave packet is studied. First 
term is f short in h and is expressed by Bessel functions. We have 

Xfshort. (107) 

Li is evaluated at a large \5t\ and we have 

Li = (7r(7.)ie'^^^* J dr e-'^^'^-^^''-'^'^''^" f short, f^xi- x^. (108) 

Here the integration is made in the space- like region A < 0. It is convenient 
to write 

ri = v,5t + fi, (109) 

and to write A in the form 

X = 5t^ - r'l - r| = 5t^ - {vjt + fif - f| ^ -2tvn5t - rf - f|. (110) 

The Li for the large \5t\ in the space-like region is written with the asymptotic 
expression of the Bessel function and becomes 



47r2 (2m) 



1 



X ( 1 ) ^irh^j2vufi\5t\+f1+f^ fill) 

V+2v^n5t + f2 + f|y ■ ^ ^ 

By the Gaussian integration around = 0, f; = —i2a^Pi,, we have the 
asymptotic expression of Li at a large \St\ 

Li = (7rcrj.)tZi, 

3 

47r2 V2m/ \+2v^2a^pJt J ' ^ ' 

Obviously Li oscillates fast as e""'^il''*l'^ where ci is determined by and a^, 
and is short-range. The integration carried out with a different stationary 
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value of ri which takes into account the last term in the right-hand side 

gives almost equivalent result. The integration of Li in the time-like region 

1 

A > is carried in a similar manner and Li decreases with time as e~'^'^i 1*^*1^ 
and final result after the time integration is almost the same as that of the 
space-like region. It is noted that the long-range term which appeared from 
the isolated 1/A singularity in Eq. (!93|) does not exist in Li in fact. The 
reason for its absence is that the Bessel function decreases much faster in 
the space- like region than 1/A and oscillates much faster than 1/A in the 
time-like region. Hence the long-range correlation is not generated from the 
Li and the light-cone singularity 6{X)e{6t) and 1/A are the only source of the 
long-distance correlation. 

Second term is from I2, Eq. f l84p . We have this term, L2, 



(27r) 



3^2, 



-i(En-Er,-^/q'2+ml-v^ip^-q-pr,))5t 



mf 



X 



^-^Ap.-'i-P.r0 (^E^ _ ^ ^2 ^ 



mf 



;ii3) 



The angular velocity of Eq. (11131) in L2 varies with q and L2 becomes to 
have a short-range correlation of the length, 2yJo\,, in the time direction. So 
the L2's contribution to the total probability comes from the small 5t region 
and this corresponds to the short-range component. 

Thus the coordinate integration of A^^^(5t, 5x) is written in the form 



2i—p^-p^ 



d 

1 + 2p^-p^—^ + 



{J5(X) + Li) + L2 



2i{TTau)^^^p^-p^ 



-'TV 



d 

l + 2p^-{-p^)—- 



3 



-i)L') 



(114) 



2 \0t\ J \ TT 

In the above equation, pi = ml is neglected since this is extremely small 
compared to m^, p^-Pu and a^. This is neglected also in most other places 
except the slow phase (f){6t). The first term in the right-hand side of Eq. (11141) 
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has the long-distance correlation and the second term has a short distance 
correlation. They are separated in a clear manner. 

4.3.3 Convergence condition 

At the end of this section, we study when our method is valid. In our 
calculations the integration region is split into the one of finite region and 
the other of infinite region and the separation of the correlation function into 
the singular term and the regular term, which is equivalent to the separation 
of the correlation function into the long-distance term and the short distance 
term, is made next. The long-distance term is generated from the light-cone 
singularity, which is obtained by the the expansion of Eq. (178|) . so for our 
method to be valid, the convergence condition and the separability of the 
long-distance behavior should be satisfied. 

We study the convergence condition when the power series 

(115) 

becomes finite using the asymptotic expression of Li, Eq. f lll2p . first. The 
most dangerous term in Li, is ma . Other terms converge when this converges. 

50 we find the convergence condition from the series 

S, = Y.(-^P.-P.r^(^)\rn^t^. (116) 
The Si becomes into the form. 

Hence the series converges in the kinematical region Eq. fjHT]) . At 2pt^ -p^, = fh? 

51 becomes finite, and the value is expressed by the zeta function, 

5i = J]n"t(^)l=C('^Vm)i (118) 
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Hence in the region, Eq. f l8ip . the correlation function has the singular terms 
and the total probability has the long-range terms Js[X) and Ji/x- In the 
outside of this region, the power series diverges and our method does not 
work. / is evaluated directly and has no long range term. The / obtained 
from the finite muon momentum is equivalent to the l2- 

We study if the power series Eq. f lllSp oscillates with time rapidly 
next. For this purpose we study 

n ^ 

When 5*2 oscillates with the other terms in Li oscillate with the 

By the explicit calculations, we have 

and the function oscillates with a/I^^I in the kinematical region Eq. (181 p . So 
the separating the light-cone singular term from Ii is valid and is used for the 
evaluation of the finite-size correction of the probability at the finite distance. 



;ii9) 



5 Time-dependent probability 



Substituting Eqs. (EHD, (171]), and ffTT4D into Eq. ([70]), we have the probability 
for measuring the neutrino at a space-time position {X,^,Ty), when the pion 
momentum distribution PexpiPir) is known, in the following form 



3 1 



X / -^Pn-Pu I dtidt2 



. m 



(27r)6"" 
Ist^i^t) 2Li 2 /a^\i 



d^ 



PexpiPir) 



\ 2 



(121) 

From the pion coherence length obtained in the previous section, the coher- 
ence condition, Eq. ([2]), is satisfied and the pion Gaussian parts are regarded 
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as constant in ti and t2, 

e-2i7(^-^-+(-^-^^^)(*-^'^)+-^-(^-^'^))' ^ constant in h, 
e-2i7(^-^-+(-^-^^^)(*-^^)+^^-(^-^'^))' ^ constant in t2, (122) 

when the integration on time ti and t2 are made in a distance of our interest 
which is of order a few 100 [m]. The integration on ti, t2 will be made in the 
next section. 

When the above conditions Eq. f ll22p are satisfied, the area where the 
neutrino is produced is inside of the same pion and the neutrino waves are 
treated coherently and are capable of showing interference. In a much larger 
distance where this condition is not satisfied, two positions can not be in the 
same pion and the interference disappears. 



5.1 Integrations on times 



Integration of the probability over the time ti and t2 are carried and proba- 
bility at a finite T is obtained now. The time integral of the slowly decreasing 
term is 



dtidt2 



where g(T, Uu) is 



2 I Sin X 



1 — cos X 



X 



X = UuT, Sin X 



^ , sin t 
dt . 



(123) 



(124) 



The slope of ^'(T, u^) at T = is 



d dx d 

— ^(T,c..)|t=o = ^^ 



—2 I Sin X — 



cosx 



X 



-cu,. (125) 



x=0 



At the infinite time T = oo, g{T,Ui,) becomes g{oo,u!iy) = — vr that is can- 
celled with the short-range term of Ii of Eq. f lTHj) . So it is convenient to 
subtract the asymptotic value from g(T,uJi^) and define g{T,Uiy) 



g{T, u^) = g(T, u^) - g{oo, u^). 



(126) 



We understand that the short-range part Li cancels with g{oo, u^) and write 
the total probability with ^(T,^;,^) and the short-range term from I2. 



36 



The time integral of the short-range term, Z/2, is 
dtidt2L2{St) 



TT V TT 



2 la, 



dtidt2 [ ,^ ^ ^-i(E^-E^-^^/q^+ml-v^-(pn-q-p^))5t 

V ^ io J ^^^(fTm^ 




X e-'^-'^P^-'^-P-'^'e [e^ - ^ql + + ml 
= TGo, (127) 

where the constant Go is given in the integral 

d^q 



^^0 = 2^^ J ^J^=^f -E,-^q^ + ml~v,- (p^ - q- )) 



X 



^-MP.-g-P.^Q (^E^ - ^qf + qf + , (128) 



and is estimated numerically. Due to the rapid oscillation in 6t, L2 con- 
tributes to the probability from the small 5t region and the integrations over 
the time becomes constant early with the time interval T. This has no finite 
size correction. The regular term Li is also the same. 



5.2 Total transition probability 

Adding the slowly decreasing part and the short-range part, we have the 
final expression of the total probability. The neutrino coordinate of the 

3 r„ 

final state is integrated in Eq.f ll2ip and a factor (o"7rVr)2 emerges. This cr^ 

3 

dependence is cancelled by the (47r/(T^)2 from the normalization Eq. (16 7p and 
the final result is independent from a^^. The total transition probability is 
expressed in the form. 



/d p-jY f d Pi, 

-^Pexp{Pn) j -^{P-K-Pu)[g{T^,UJv) + Gq], 

Do = \N^.?MHo^'n)^^^ = ^ (129) 

where L = cT is the length of decay region. Eq. (11291) depends on the neutrino 
wave packet size a,. 
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Fig. 3: The neutrino probability in the forward direction per time at a finite 
distance L is given. The constant shows the short-range normal term and 
the long-range term is written on top of the normal term. The horizontal 
axis shows the distance in [m] and the probability is of arbitrary unit. Clear 
excess is seen in the distance below 1200 [m] . The neutrino mass, pion energy, 
neutrino energy are 1 [eV/c^], 4 [GeV], and 800 [MeV]. 



5.2.1 Neutrino angle integration 

In the normal term Gq of Eq. f ll29p the cosine of neutrino angle 6 is deter- 
mined approximately from 

{p.-p.y=pl = ml, (130) 

because the energy and momentum conservation is approximately satisfied 
in the normal term. Hence the product of the momenta is expressed by the 
masses 

2 2 

— m 

Pn-Pu = — (131) 
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Fig. 4: The total probability integrated in the neutrino angle per time at 
a finite distance L is given. The constant shows the short-range normal 
term and the long-range term is written on top of the normal term. The 
horizontal axis shows the distance in [m] and the probability is of arbitrary 
unit. The excess becomes less clear than the forward direction, but is seen 
in the distance below 1200 [m]. The neutrino mass, pion energy, neutrino 
energy are 1.0 [eV/c^], 4 [GeV] and 4.5 [GeV] , and 800 [MeV]. 



and the cosine of the angle satisfies 



1 — cos 9 



ml - ml 



ml 



m\\pu\ m?- ^^^^^ 

The cos^ is very close to 1. On the other hand, the long-range component of 
the neutrino probability, g{T,Uiy) of Eq. (11291) . is derived from the light-cone 
singular term. This term is present only when the product of the momenta 
is in the convergence domain Eq. (1811) . Hence the long-range term is present 
in the kinematical region. 



m: 



m. 



(133) 
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Since the angular region of Eq. 



is slightly different from Eq. ( 11311) 



and it is impossible to distinguish the latter from the former region experi- 
mentally, the neutrino angle is integrated. We integrate the neutrino angle 
of both terms separately. We have the normal term, Gq, in the form 



E, 




Tc \a 



dq 



m, 



X 



{2n) 



2 I ml -ml 



El/, max 



where 



E,. 



ml - ml 



E,, 



dE„ 



mz 



m. 



(134) 



(135) 



2{E^ + \p^\y ^'"^"^^ 2(K-|p-;i) 

and the Gaussian function is approximated by the delta function for the 
computational convenience. The angle is determined uniquely. We have for 
the long-range term, g{T,u„), in the form 



E, 



2n 



2tt 



2n 



\Pu\^d\p^ 
E, 

\Pu?d\p^ 



d COS 6{E^E^ - \Pj,\\Pu\ cos6)g{T,Ui, 



E,/ 



El/, max 



dE,/ 



E-^Ey COS 6^ - i Iptt I \Pv I cos^ 6 



\Vn\\Vv\ 



2\p. 



iE.E^ - mvAf ^(T,c^.), (136) 



where the angle is very close to the former value but is not uniquely deter- 
mined. Finally we have the energy dependent probability 

dP 



dE„ 



d^Pn / - N 27r 

Pexp\P-K) ■ - ■ X 



E. 



\pA 



m, 



(137) 
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5.3 Comparison with measurements 
5.3.1 Sharp pion momentum 

When initial pion has a discrete momentum P^r, the PexpiPir) is given as 

PexpiPir) =5{Pn- -Pyr) (138) 

and the energy- dependent probability for this pion is 



dP 2 2 n 1 27r 



7r(m^ - m^ 



(139) 



Eq. f ll39p is independent from the position X^, of the initial state and depends 
upon the momenta and the time difference T = T,y — T^^. Consequently 
although quantum mechanics tells that the probability Eq. fll39p should be 
applied to events of one initial state of the pion, it is possible to average over 
different positions. The result is obviously the same Eq. fll39p . 

At T — > oo, g(T,oj^) vanishes and the decay rate per time of the energy 
i?7r is given from the first term of Eq. f ll39p as 

1 27]- rEv,max 

P/T = g'mlD, — -^Mml - ml) / dE, 

The stationary value is independent from the wave packet size and is consis- 
tent with [21]. This value furthermore agrees to the standard value obtained 
using the plane waves. Hence the large time limit of our result is equivalent 
to the known result. 

At finite T, the diffraction term does not vanish and Eq. fll39p deviates 
from the standard value obtained using the plane waves. The relative mag- 
nitude of the diffraction term ^(T, Uv) to the normal term Gq is independent 
from detection process. So we compute g{T,Uy) and Go of Eq. (11290 at the 
forward direction = and the energy dependent total probability that is 
integrated over the neutrino angle in the following. The function ^(T, u) and 
Go are plotted in Fig. [3] for the mass of neutrino, m,^ = 1 [eV/c^], the pion 
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Fig. 5: The total probability integrated in the neutrino angle per time at a 
finite distance L is given. The constant shows the short-range normal term 
and the long-range term is written on top of the normal term. The horizontal 
axis shows the distance in [m] and the probability is of arbitrary unit. Clear 
uniform excess is seen in the distance below 1200 [m]. The neutrino mass, 
pion energy, neutrino energy are 0.6 [eV/c^], 4 [GeV], and 800 [MeV]. 

energy i?^ = 4 [GeV] , and the neutrino energy E^^ = 800 [MeV] . For the wave 
packet size of the neutrino, the size of the nucleus of the mass number A, 
ay = A3 /ml is used. The value becomes = 6.4/m^ for the ^^O nucleus 
and this is used for the following evaluations. From this figure it is seen that 
there is an excess of the fiux at short distance region L < 600 [m] and the 
maximal excess is about 0.4 at L = 0. The slope at the origin L = is 
determined by u^. The slowly decreasing term that is generated from the 
singularity at the light cone has a finite magnitude. 

The total probability that is integrated over the neutrino angle Eq. fll37p 
is presented next. The probability for the neutrino mass my = 1.0 [eV/c^] 
and the pion energy 4 [GeV] and 4.5 [GeV] are given in Fig. HI and for the 
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Fig. 6: The total probability integrated in the angle per time at a finite 
distance L is given. The constant shows the short-range normal term and 
the long-range term is written on top of the normal term. The horizontal axis 
shows the distance in [m] and the probability is of arbitrary unit. Clear excess 
and decreasing behavior are seen in the distance below 600 [m]. The neutrino 
mass, pion energy, neutrino energy are 1 [eV/c^], 4 [GeV], and 100 [MeV]. 



smaller neutrino mass m,^ = 0.6 [eV/c^] is given in Fig. [51 Gq is unchanged 
with the distance but the long-distance term, ^(T,^;,^), decreases slowly with 
the distance than that of rrii, = 1 [eV/c^]. Hence the longer distance is nec- 
essary if the mass of the neutrino is even smaller. For the muon neutrino, 
it is impossible to measure the event at a energy lower than few 100 [MeV]. 
The electron neutrino is used then. Considering the situation for the electron 
neutrino, we present the total probability for the lower energies. The prob- 
ability for the neutrino mass mj, = 1.0 [eV/c^] with the energy 100 [MeV] is 
given in Fig.O The slowly decreasing component of the probability becomes 
more prominent with lower values. Hence to observe this component, the 
experiment of the lower neutrino energy is more convenient. 



43 




400 600 800 

Length of decay volume [m] 



1200 



Fig. 7: The neutrino probability integrated in the neutrino angle per time 
at a finite distance L is given. The constant shows the short-range normal 
term and the long-range term is written on top of the normal term. The 
horizontal axis shows the distance in [m] and the probability is of arbitrary 
unit. Clear excess is seen in the distance below 1200 [m]. The neutrino mass, 
pion energy, neutrino energy are 0.1 [eV/c^], 4 [GeV], and 10 [MeV]. 



The typical length Lq of the diffraction term is 



Lo [m] = ^ = m -:^'Z:L - (141) 



By the observation of this component together with the neutrino's energy, the 
determination of the neutrino mass may becomes possible. The neutrino's 
energy is measured with uncertainty /S.Ey, which is of the order 0.1 x Ey. This 
uncertainty is 100 [MeV] for the energy 1 [GeV] and is accidentally same order 
as that of the minimum uncertainty h/5x derived from Eq. ( I20l) . The total 
probability for a larger value of energy uncertainty is easily computed using 
Eq. f ll29p . Figs. (3)-(7) show the length dependence of the probability. If the 
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Fig. 8: The energy dependence of the probabihty integrated in the angle 
at distance L = 100 [m] is given. The lower curve shows the short-range 
normal term and the long-range term is added on top of the normal term. 
The horizontal axis shows the neutrino energy in [MeV] and the probability 
is of arbitrary unit. The neutrino mass and pion energy are 1.0 [eV/c^] and 
4 [GeV]. 

mass is around 1 [eV/c^] the excess of the neutrino flux of about 20 percent 
at the distance less than a few hundred meters is found. In the long-baseline 
neutrino oscillation experiments, the neutrino flux at the near detectors has 
observed excesses of about 10 — 20 percent [281 EHl [30] . We believe this is con- 
nected with the excesses found in this paper. We use mainly = 1 [eV/c^] 
throughout this paper. Because the probability has a constant term and the 
T-dependent term, the T-dependent term is extracted easily by subtracting 
the constant term from the total probability. The slowly decreasing com- 
ponent decreases with the scale determined by the neutrino's mass and the 
energy. Although the excess of the flux would be found, the decreasing be- 
havior becomes difficult to observe if the mass is less than 0.1 [eV/c^] for the 
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muon neutrino. In this case, the electron neutrino is useful. The electron 
neutrino is produced in the decays of the muon, neutron, K-meson, and nu- 
cleus. In these decays the present mechanism works. So we plot the figure for 
rriu = 0.1 [eV/c^], Ei, = 10 [MeV] in Fig. [3 A decreasing part is clearly seen. 
So in order to observe the slow decreasing behavior for the small neutrino 
mass less than or about the same as 0.1 [eV/c^], the electron neutrino should 
be used. The decay of the muon and others will be studied in a forthcoming 
paper. In Fig. |8] the energy dependence of the total probability is given. The 
long-range term has a different form and has a maximum at the neutrino 
energy ^ l/3E^^max- 

5.3.2 Wide distribution of pion momentum 

When the momentum distribution Pexpiv-K) of the initial pion is known, the 
energy-dependent probability of the neutrino is given by Eq. f ll37p . Eq. f ll37p 
is also independent from the position X^^ and depends upon the momentum 
of the pion and of the neutrino and the time difference T = T,^ — T^. In 
the experiments, the pion's position is not measured and the average over 
the position is made. This average probability agrees to Eq.f ll37p . The 
probability varies slowly with the pion's momentum and is regarded constant 
in the energy range of order 100 [MeV]. So the experimental observation of 
the diffraction term is quite easy. 

5.4 On the universality of the long-distance correction 

The long-distance correction of the neutrino detection probability was ob- 
tained using the wave packets. The correction is decreasing extremely slowly 
with the distance. The slope is determined only from the mass and energy 
of the neutrino and is independent from details of other parameters of the 
system such as the size or shape of the wave packets and the position of 
the final wave packets and others. Hence this is the genuine property of the 
wave function |muon, neutrino (t)), of Eq. fl42p . The slope of this term has 
the universal property and is capable of measuring in experiments. 

5.4.1 Muon in pion decay 

When the muon is observed in the same processes, the anomalous behavior 
is determined by the mass and energy of the muon as w?^/2E^. Since the 
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muon mass is larger than the neutrino mass by 10^, the typical length for the 
muon is smaller than that of the neutrino by 10~^^. For the muon of energy 
1 [GeV], the length is order 10~^^ [m]. This value is a microscopic size. So the 
muon probability becomes constant in this distance and has no observable 
long-range effect. The muon from the pion decay has no coherence effect and 
can be treated incoherently. 

6 Summary and implications 

In this paper, we studied the interference of the neutrino in high energy 
pion decay. The position-dependent neutrino flux was computed and the 
new long-distance component of universal property that is insensitive to the 
pion's initial conditions was found to exist in addition to the normal in- 
coherent term. The pion's coherence length estimated in Section 2 shows 
that this is long enough for the new universal term to be observed. Since 
the new term reflects the interference of the neutrino, it decreases with the 
distance in the universal manner that depends upon the absolute value of 
the neutrino mass. 

The present interference phenomenon is caused by the light-cone singu- 
larity of the correlation function of the pion and muon system A^_^((5t, (5a;) 
and the slow phase 0c of the neutrino wave function, Eq. flSHl) These two 
ingredients are characteristic phenomena of the relativistic invariant system. 
The light-cone singularity is formed from the waves of the infinite momen- 
tum, which have always the light velocity. The slow phase of the neutrino 
wave packet is the outcome of the cancellation between the time and space 
oscillations. The phase is determined by the difference of space-time co- 
ordinates {St, Sx) and the central values of the energy and momentum, as 
E(p)6t — p-6x, where the energy -E'(p) is given by a/p^ + m^. When the 
difference of positions 6x is moving with the light velocity in the parallel 
direction to the momentum p, = c{p/p), as 6x = VcSt, then the total phase 
becomes Ey5t—pu-5x = {m1/2Ey)5t. The angular velocity becomes the small 
value ml/2Ei, due to the relativistic invariance and makes the interference 
phenomenon long-range. The new term in the probability decreases slowly 
with the distance in the universal manner determined by the mass and en- 
ergy of the neutrino as m1/2Ey. This form is independent from the details of 
wave packet shape and parameters as far as the reality of the neutrino wave 
function is satisfied, which is ensured from the invariance under the 
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time inversion. The relative magnitude of this component is not universal 
and depends upon the size of wave packet. Based on the estimation of the 
size, we found that the magnitude of the new universal term is sizable for the 
measurement. Since the slope of the probability, Eq. fll29p . is determined by 
the mass and energy of the neutrino, the absolute value of the neutrino mass 
would be found from the neutrino interference experiments. 

The new component shall be observed as the excesses of the neutrino 
flux in the ground experiments. The excesses of the neutrino flux at the 
macroscopic short distance region of the order of a few hundred meters were 
computed and shown in Figs. (3)- (8). From these figures, the excesses are 
not large but are sizable magnitudes. Hence these excesses shall be observed 
in these distances. Actually fluxes measured in the near detectors of the long- 
baseline experiments of K2K ^28j and MiniBooNE [29] may show excesses of 
about 10 — 20 percent of the Monte Carlo estimations. Monte Carlo estima- 
tions of the fluxes are obtained using naive decay probabilities and do not 
have the coherence effects we presented in the present work. So the excess of 
these experiments may be related with the excesses due to interferences. The 
excess is not clear in MINOS j30j. With more statistics, qualitative analysis 
might become possible to test the new universal term on the neutrino flux 
at the finite distance. From Figs. (3)-(8), if the mass is in the range from 
0.1 [eV/c2] to 2 [eV/c^], the near detectors at T2K, MiniBooNE, MINOS and 
other experiments might be able to measure these signatures. The absolute 
value of the mass could be found then. It is worthwhile at the end to sum- 
marize the reasons why the interference term of the long-distance behavior 
emerges in the pion decay and is computed with the wave packet represen- 
tation. The connection of the long-distance interference phenomenon of the 
neutrino with the Heisenberg's uncertainty relation is also addressed here. 

Due to relativistic invariance, the correlation function A^_^((5t, 6x) has a 
singularity near A = 0, which is extended to large distance \Sx\ — )• oo, as 
was shown in Eq. (!82|) . This is one of the features of relativistic quantum 
fields and is one reason why the long-range correlation emerged. For a non- 
relativistic system, on the other hand, the same calculation for stationary 
states is made by. 



and the only one point 6x = satisfies the condition. Long-range correlation 
is not generated. The rotational invariant three-dimensional space is com- 




(142) 
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pact but the Lorentz invariant four-dimensional space is non-compact. So 
it is quite natural for the non-relativistic system not to have the long-range 
correlation that the relativistic system has. The light-cone singularity is the 
peculiar property of the relativistic system. 

Heisenberg uncertainty relation is slightly modified in the wave along the 
light cone. The neutrino wave function behaves at large-distance along the 
light-cone region in the form 

m? 

Mt, x) = / = f—-, (143) 

X ct 

where / has no dependence on the distance \x\. Consequently the uncertainty 
relation between the energy width 6E and the time interval 6t becomes 

5t5^ = St5E x-(^X ^h. (144) 

The ratio {m^/Eyf is of the order 10~^^ and 5t becomes macroscopic even 
if the energy width 5E is microscopic of order 100 [MeV]. For instance if the 
pion Compton wave length, A^, is used for the microscopic length, then c5t 
becomes 

10^^ lO^m, (145) 

which is about the distance between the pion source and the near detector 
in fact. So interference effect of the present paper appears in this distance 
and is observable using the apparatus of much smaller size. 

The position-dependent probability was computed easily with the use of 
the wave packet representation. Any representation is equivalent to plane 
waves as far as the whole complete set of functions is used. The ordinary 
probabihty of the transition process is defined from the states t — ±oo. 
The normal scattering amplitude is the overlap between the in-state aX, t — 
— oo and out-state ai t — oo, and the space-time coordinates are integrated 
from — cxD to oo and the energy and momentum of the final state is the same 
as that of the initial state. Hence in the pion decay, the momentum of the 
muon in the final state of the ordinary scattering experiments are bounded. 
So the infinite momentum is not included in the muon of the final state. 
From these amplitudes, the amplitudes and probability at the finite-time 
interval are neither computable nor obtained. In the wave packet formalism. 
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on the other hand, it is possible to compute the amphtude and probabihty at 
finite-time interval directly. Energy and momentum conservation is slightly 
violated in this amplitude and the states of infinite momentum couple and 
give the light-cone singularity to the correlation function A^^^(5a;). The 
contribution from these states vanishes at infinite time and these states do 
not contribute to the probability measured at infinite distance, i.e., ordinary 
S-matrix. Thus the light-cone singularity gave the interference term, which 
ended up as the finite-time interval effect of the neutrino probability. Wave 
packet formalism gives new information that is not calculable in the standard 
scattering amplitude. Hence our calculation does not contradict with the 
ordinary calculation of the S-matrix in momentum representation but has 
advantage of giving a new universal physical quantity directly. 

Since characteristic small phase of the relativistic wave along the light 
cone and the light-cone singularity are derived from Lorentz invariance, there 
would be similar effect if there exists other light particle. Axion is a possible 
candidate of light particle and might show the long-distance interference 
phenomenon. 

In this paper we ignored higher order effects such as pion life time and 
pion mean free path in studying the quantum effects. We will study these 
problems and other large scale physical phenomena of low energy neutrinos 
in subsequent papers. 
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Appendix A Light-cone singularity 

A-I Light-cone singularity and small neutrino mass 

A conjugate momentum to {6t, 6x) is p,r ~P^l from Eq. f lTTj) and the invariant 
square of this momentum becomes 

(Ptt - P^Y = ml + ml - 2 (^^/ff+m^^ {^^Jpl + mj^ + 2\p^\\p ^\ cos 9 , 

(146) 
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where 6 is the angle between the pion and muon momenta. The invariant 
vanishes when the cosine of angle becomes 




+ small terms. 



(147) 

This equation has a solution in a finite muon momentum region for a given 
pion momentum. 

A-II Long-range correlation for general wave packets 
A-II.l Non-Gaussian wave packet 

Non-Gaussian wave packets were studied in the general manner in the text 
and the universal behavior of the phase was obtained. In this appendix the 
explicit forms of the wave packets are studied. It is reconfirmed that the 
long-range component of the probability at around t = 2TTEy/m'^ becomes 
the universal form, 
type 1 

One way to express the non- Gaussian wave packet is to multiply Hermi- 
tian polynomials and to write the amplitude in the form 

j dkye-'^^^'^-^''^"H^{^,{K-p,))e'^''^^^^^^ (148) 
(27r)2 J 

where Hn is assumed to be real in order the wave packets to preserve the 
time reversal symmetry and an even function of kj, — p^, in order the wave 
packets to preserve parity , Eqs. ( l28l) and ( !32l) . 

Since we study symmetric wave packets, it is sufficient to prove the sim- 
plest case 

Hr, = a,{K-p,f. (149) 

The spreading effect was studied in the previous appendix and does not 
change the final result. So we ignore the spreading effect here. The momen- 
tum integration Eq. ( |49|) is replaced with 
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After the straightforward calculations we have this integral in the form, 

ku-Pu-\ {x - Xy- Vy{t -Ty)}\ -{x - Xy - Vu{t - Ty)Y 



Of. 



+ — [K-PU + —{X- Xy - Vu{t - Ty)} ] {X - Xy- Vu{t - Ty)} 



I I ^i(E{p^)(t~T„)-p^-{x-X^))-,^(x-X^-v^{t-T„)f 



X 3 {X-Xy-Vy{t-Ty)}' 



(151) 



Next the integral Eq. (IHHjl is studied. This becomes for the non-Gaussian 
wave packet to the integral 



Jsix)=K / dx^dx2e'^'''^>'^ 



i^{5x)-^(xi-X^~v,,{ti-n)f-^{x2-X^-v^{t2-n)f 



5(A) 



X 3 {fl -Xy- Vy{h - Ty)y 3 { X 2 " X, - ^7,(^2 " Ty)} 



O-y 



Oy 



and is written by using the center coordinate X^ = 
dinate = x^ — X2 in the form 



(152) 

and relative coor- 



A^' / dXdfe 



,^^Sx)-^X^-^r-^' 



li. 



1 



9 _ _ 2X^ + -f M + — - -X'f' + —f 
cr\ 2 al\ 2 16 



1 



(153) 



where 



X = X-Xy- Vyi^X"^ - Ty), F= T " VT^. 

The integration on X and r are made and we have the final result 

Jsw = N'^Mlay^^e^^^-^^>° 

-f + ^^(i--.)^(r°)^ + o((i-..)VT) 



(154) 
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Thus the phase factor has the same universal form as the Gaussian wave 
packet and the correction is determined by the small parameter {l—v^)'^{r^y 
in the form 



.2 \ 2 



hence the correction is negligible at high energy. 

We have proved that the correlation function of the non-Gaussian wave 
packet has the same slow phase and long-range term as the Gaussian wave 
packet and the small correction becomes negligible for the simplest case 
Eq. ll49[ Hence for any polynomials Hn that are invariant under the time 
and space inversions, the correlation function has the same long-range term 
and small negligible corrections. 

type 2 

Another way to write the non-Gaussian wave packet is to use a function 
a{p), and to write 



2 



^ / t/4e-"(^^'')+<^(^^-)(*-^'')-^^-("^-^-)). (156) 



The large t = T behavior is found by the stationary momentum which satis- 
fies the equation 

Symmetric real wave packet is assumed also here from parity and time re- 
versal invariances of the wave packets and we write, 

a{K) = a(p.) + ^^^^^^a. + {K - Pu)-{K - Pu)ptJ + ■ ■ ■ , (158) 

where the a and Cij are real numbers. The momentum integration of Eq. (11561) 
becomes the form 

(27r)i J " 
(27r)i J ' 

x(l - {{k, - p^)^)\{k, - p,),fCij). (159) 
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The correction to the Gaussian wave packet is generated by the higher order 
terms of k^, — p^, in the right hand side and is treated in a same way as the 
previous type 1 case. Hence this integral has the leading long-range term 
which is equivalent to that of the Gaussian wave packet and the negligible 
correction expressed by Eq. f llSSp . 

For studying the asymptotic behavior at t — T^, — )■ oo we solve the sta- 
tionarity equation, 



d 



a{k,)-i{E{k,){t-T,)-k,-{x-X,)} =0 (160) 



and expand the integral around the stationary momentum. The wave in the 
transverse direction to this momentum spreads but spreading is very small 
in the longitudinal direction [TU]. From the result of the previous appendix, 
the final result is the same and so is not presented here, 
type 3 

In the type 1 and 2 the time reversal and parity symmetries are assumed 
for the wave packet shape. If these symmetries are not required, the function 
Hn or a has an imaginary part. In this case, the correlation function has a 
correction term in the order (1 — f )(ti — t2) and this term is expressed 

1 TTL^ 

{l-v,){h-h) = —-^{h-t,) (161) 

hence the correction term vanishes at the high energy. With a suitable param- 
eter, the universal form of the slowly decreasing component of the probability 
of the present work may become observable even in arbitrary system. The 
Lorentz invariant form of the energy dependent phase of the wave packet 
and the light-cone singularity of the pion and muon decay vertex give this 
universal behavior. 
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